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Abstract 

Y-Al 2 0 3 is a porous metal oxide and described as a defective spinel with some cationic vacancies. In this work, we 
calculate the electronic density of states and band structure for the bulk of this material. The calculations are 
performed within the density functional theory using the full potential augmented plan waves plus local orbital 
method, as embodied in the WIEN2k code. We show that the modified Becke-Johnson exchange potential, as a 
semi-local method, can predict the bandgap in better agreement with the experiment even compared to the 
accurate but much more expensive green function method. Moreover, our electronic structure analysis indicates 
that the character of the valence band maximum mainly originates from the p orbital of those oxygen atoms that 
are close to the vacancy. The charge density results show that the polarization of the oxygen electron cloud is 
directed toward aluminum cations, which cause Al and 0 atoms to be tightly connected by a strong dipole bond. 
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Background 

y-Al 2 0 3 is an important material in nanotechnology be- 
cause of its porous structure, high surface area, and high 
catalytic surface activity. This material is widely used as 
a catalyst, an adsorbent, and a support for industrial cat- 
alysts in the oxidation of organics and the catalytic re- 
duction of automotive pollutants such as nitric oxide [1]. 
This oxide material is a metastable 'transition' phase of 
aluminum oxide or alumina which can be formed via de- 
hydrating boehmite y-AlOOH at rather low tempera- 
tures (350°C to 1,000°C) [2] and is able to keep its 
crystalline structure unchanged up to about 1,200°C [3]. 
Recently, this material is considered as a suitable alter- 
native to silicon for producing semiconductor nonvola- 
tile random access memories for future applications [4]. 
Thus, investigation of the atomic and electronic struc- 
tures of y-Al 2 0 3 has also received much attention over 
the last decade [5,6]. 

The density functional theory (DFT)-based methods 
are appropriate to investigate the electronic structure of 
solid materials. However, the precise electronic bandgap 
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prediction is a shortcoming in the DFT-based methods 
(see Table 1). 

In this work, we have calculated the electronic struc- 
ture of y-Al 2 0 3 compound using the mBJ exchange po- 
tential [7]. Our result shows that the semi-local mBJ 
method predicts the y-Al 2 0 3 bandgap better than the 
LDA and GGA when compared to the experimental 
data. The results which are obtained by the semi-local 
mBJ method are in good agreement with those of the ex- 
perimental data. In fact, the precision of theoretical pre- 
dictions in the mBJ scheme is comparable with the 
accurate but much more expensive Green function 
(GW) method. In the mBJ exchange potential, there is a 
correction (c)-factor which can be used as an adjustable 
parameter for improving the bandgap prediction. There- 
fore, the mBJ method has the capability to overcome the 
well-known shortcoming of the DFT-based methods in 
predicting the bandgaps. In this paper, we have utilized 
this capability to yield the bandgap of y-Al 2 0 3 by deter- 
mining and fixing the adjustable c parameter of the mBJ 
potential. The calculated bandgap within our fixed mBJ 
calculations is found to be in better agreement with that 
of the experiment than that of the GW method. 
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Table 1 Bandgaps 



References 


Method 


Bandgap (eV) 


LoJ 


r r a 


4.U 


[14] 


LDA 


4.2 




PW91-GGA 


3.8 


[13] 


G 0 Wo 


7.2 


Current 


GGA 


4.13 




Regular mBJ 


6.19 




Non-regular mBJ 


8.02 


Exp. [12] 




8.7 



GGA, generalized gradient approximation; LDA, local density approximation; 
GW, Green function; mBJ, modified Becke-Johnson. 



Main text 

Methods 

All of the calculations in this work were carried out 
using the WIEN2K code [8], which is based on the full 
potential augmented plane waves plus local orbital 
method. For the calculations, 96 irreducible points in 
the first Brillouin zone were used. The points were sepa- 
rated for integration over the Brillouin zone of the unit 
cell using a set of 8 x 8 x 3 special k-points. The expan- 
sion of the wave functions and charge densities was cut 
off by setting the RmtKmexj G Max , and Muffin-tin radii, 
as tabulated in Table 2. In this research, we have used 
the mBJ exchange potential to investigate the electronic 
structure and bandgap of bulk y-Al 2 0 3 . This functional 
has been proposed by Tran and Blaha [7]. It is a modi- 
fied version of the BJ exchange correlation potential 
(introduced by Becke and Johnson [9]) to improve the 
bandgap prediction of the DFT-based methods. In the 
regular mBJ calculations, we are only allowed to deter- 
mine the c-factor self-consistently through our ab initio 
calculations. It is known that the mBJ potential cannot 
always successfully predict the bandgaps for all of the 
materials in the world in excellent agreement with the 
experiment [7]. However, it may be possible to efficiently 
improve the bandgap perdition by performing non- 
regular mBJ calculations. The c-factor is self-consistently 
determined in the regular mBJ calculations. In the non- 
regular mBJ calculation, one can internally increase the c 
parameter of the mBJ method to obtain better results. It 
has been shown that the bandgaps increase as the c par- 
ameter increases [7,10]. Thus, one can perform a non- 
regular mBJ calculation by increasing the c parameter to 
an optimized c value. In order to more precisely deter- 
mine the bandgap for the y-Al 2 0 3 compound, we 
optimize the c-factor and then adjust the c-factor 



Table 2 Computational parameters 



K-point 
meshes 


Rmt-Kmax 
(bohrxRy 172 ) 


Separation 
energy (Ry) 


Giviax 

(Ry 1/2 ) 


Muffin-tin 
radii (bohr) 


8x8x3 


7 


8 


12 


1.51 



manually to the optimized value in our non-regular mBJ 
calculations (see [10]). We show that the non-regular 
calculation can efficiently improve the bandgap. 

Crystal structure 

The perfect crystalline structure of y-Al 2 0 3 is described 
as a defective spinel, denoted as □ 2 2A/ 21 i0 3 2, (□ = Al va- 
cancy) [5,11]. To understand this formula, we should 
consider a perfect spinel, such as the MgAl 2 0 4 structure. 
However, due to stoichiometry, y-Al 2 0 3 is not a perfect 
spinel because its cation/anion ratio is less than a stoi- 
chiometric spinel structure (i.e., 24/32). Hence, 2| 
aluminum vacancies are formed in each cell. Aluminum 
vacancies can present both in tetrahedral and octahedral 
sites, as the cations are divided in tetrahedral and octa- 
hedral positions in a stoichiometric spinel. However, the 
lowest energy configuration of y-Al 2 0 3 occurs when all 
of the aluminum vacancies are in octahedral sites with 
the largest possible inter-distances. The space symmetry 
group of a perfect spinel, such as the MgAl 2 0 4 structure, 
is Fd3m. In this symmetry group, the primitive cell is a 
triclinic cell of 14 atoms, and hence, it is denoted as 
Mg 2 Al 4 0 8 . An ideal spinel unit is therefore constructed 
with the lattice constant 7.911 A [5], and all of the Mg 
atoms are replaced by aluminum. The cubic unit cell 
contains three Al 6 O g primitive cells on 'top' of each 
other, and as a result, the unit cell contains 18 Al and 24 
O atoms. To get the correct stoichiometry, we subse- 
quently remove two Al atoms from the octahedral sites 
which have the largest inter-distances. In this way, a 40- 
atom triclinic cell containing eight A1 2 0 3 formula units 
and two Al vacancies is obtained. In order to prepare a 
hexagonal cell similar to the structure already reported 
by Pinto and co-workers [6], we change the basis 
vectors as a' = b,b' = b. c' = a + b — c, where a = 
7.91l(0ii)> = 7.91l(i0i),c = 7.911(||0) . 

The prepared cell is then used for performing calcula- 
tions. For a better illustration of the Al vacancy posi- 
tions, a y-Al 2 0 3 supercell is constructed by doubling the 
size of the cell, as shown in Figure 1. 

Discussion 

Density of states 

Total density of states (DOSs) is shown in Figure 2 
within the mBJ (Figure 2a) and GGA (Figure 2b) for the 
defective spinel structure of y-Al 2 0 3 . The total DOSs 
show that the mBJ causes the unoccupied states to move 
away from the Fermi level towards higher energies with 
respect to the conduction states calculated by the GGA. 
As can be seen from Figure 2a,b, the calculated occupied 
states within both the mBJ and GGA touch but do not 
cross the Fermi level. Therefore, according to the mBJ 
prediction, the valence and conduction states repel each 
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other more strongly than that of the GGA. Hence, the 
stronger mBJ repulsion potential gives a larger bandgap 
than that of the GGA in more agreement with that of 
the experiment [12]. However, the predicted bandgap by 
the mBJ still underestimates the experimental value. We 
will discuss how to improve the bandgap prediction of 
the mBJ method in the subsequent section. 

Partial p-orbital DoSs (partial p-DOSs) are calculated 
for each of the nonequivalent oxygen and aluminum 
atoms using mBJ in the y-Al 2 0 3 compound. The p- 
DOSs of the nonequivalent O atoms are added to each 
other to obtain total p-O-DOS, as shown in Figure 2c. 
Similarly, total p-DOS for the Al atoms is the sum over 
the partial p-DOSs of the nonequivalent Al atoms 
(Figure 2d). The total p-DOSs for O and Al atoms show 
that contribution of the oxygen p-states dominates 
valence states in the vicinity of the Fermi level. There- 
fore, most of the occupied states which are touching the 
Fermi level are constituted by the p-O-DOS. This im- 
plies that the oxygen atoms and, in particular, their 
p-sates are more important for tuning the band gap 
of y-Al 2 0 3 . The latter result is deduced from the total 
p-DOS, obtained by adding p-DOSs of the nonequivalent 
O atoms. Hence, it is not still clear which of the oxygen 
atoms play a more important role. To clarify which of 
the oxygen atoms in the y-Al 2 0 3 compound are more 
appropriate for engineering the bandgap, we perform 
more elaborations on the electronic structure of the sys- 
tem by classifying the O atoms in two different groups. 
The first group contains the O atoms that are far from 
the vacancies. In the second group, the O atoms are 
selected to be close to the vacancies. The partial p-DOSs 
of the first and second groups are shown in Figure 2e,f, 
respectively. The results elucidate that the oxygen atoms 
in the second group play a more important role in the 
engineering of the bandgap, as a contribution of the 



valence p-DOS of those of oxygen atoms which are 
closer to the vacancies predominates in touching the 
Fermi level. 

Electronic structure and bandgap 

Our calculation using the regular mBJ method yields a 
bandgap of 6.195 eV, as shown in Table 1. Although this 
value is much better than the previously calculated 
results obtained by the local and local-like density 
approximations, i.e., LDA or different versions of GGA, 
as can be seen from Table 1, it still underestimates the 
experimental value compared to that of the Green func- 
tion method. In order to improve the bandgap predic- 
tion of the mBJ method, we calculate the band structure 
within a non-regular mBJ method (Figure 3b) in addition 
to the GGA band structure (Figure 3a) and the regular 
mBJ band structure. As can be seen from Figure 3b, the 
bandgap is significantly improved by our non-regular 
mBJ calculation compared to the GGA calculation 
(Figure 3a) and regular mBJ calculation (not shown), 
when compared to the experimental value of 8.7 eV [12]. 
More interestingly, our non-regular mBJ result (8.02 eV) 
is much better than that of the accurate G 0 W 0 method 
(7.2 eV) [13], compared to the experimental value 
(8.7 eV) [12]. It is worth noting that the cost of GW cal- 
culation is more expensive than that of the mBJ since 
the mBJ is still a semi-local approximation. We now dis- 
cuss the non-regular mBJ calculation. The regular mBJ 
prediction can be improved by gradually increasing the c 
parameter internally from its self-consistently converged 
value. The c parameter is converged to 1.45 for y-Al 2 0 3 
within the regular mBJ calculation. However, this value 
is so small to be used for reproducing the experimental 
bandgap. Therefore, we increase the c parameter step by 
step. In this way, we optimize the c-factor of the mBJ 
method for our case. The optimized c-factor is found to 
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Figure 2 Total and partial DOSs for y-AI 2 0 3 . Total DOSs for y-AI 2 0 3 by (a) mBJ and (b) GGA. Total partial p-DOSs within mBJ by adding all the 
p-DOSs of the nonequivalent (c) 0 atoms and (d) Al atoms, in the supercell shown in Figure 1b. Partial p-DOSs for those of 0 atoms which are 
(e) far from and (f) close to the vacancy. Fermi level is set to zero. 



be 1.8 for the y-Al 2 0 3 compound. In our non-regular 
mBJ calculations, we use this optimized value for the c 
parameter, keeping the parameter fixed during the whole 
of the self-consistent iterations, as discussed in [10]. This 
type of c-fixed mBJ calculation gives a direct bandgap of 
8.02 eV at the T point, which is closer to the experimen- 
tally measured value of 8.7 eV [12] compared with those 
of both the regular mBJ and GW methods. 

The bandgap of the system is improved by performing 
the non-regular mBJ calculation. However, it is 
mandatory to examine whether the internal structure of 
bands is destroyed by increasing the strength of the re- 
pulsion mBJ potential. To this end, one needs to com- 
pare Figure 3a and Figure 3b. The comparison 
authenticates that the band structure is not internally 



affected by the stronger repulsion potential. Indeed, the 
conduction bands are altogether shifted up by the 
same value after increasing the c parameter, keeping 
their previous structures almost the same as before. This 
is consistent with our discussion presented earlier con- 
cerning the behaviors of the regular mBJ- and GGA- 
DOSs, as can be rechecked by comparing Figure 2a and 
Figure 2b. 

Physical backbone of increasing the c-factor for the 
defective spinel y-AI 2 0 3 

In order to physically support our success in predicting 
the bandgap within the non-regular mBJ calculations in 
better agreement with that of the experiment than that 
of the GW prediction, as presented in the preceding 
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subsection, we attempt to discuss the physical mechan- 
ism underlying such an agreement in this subsection. In 
fact, the bandgap is well reproduced for the defective 
spinel y-Al 2 0 3 because the c-factor is increased to an 
optimized value. Here, we discuss that such increasing 
and optimization are necessary for the defective spinel 
y-Al 2 0 3 compound due to its existed aluminum vacan- 
cies, as aforementioned in the 'Crystal structure' subsec- 
tion. To do this, it is useful to make clear the role of the 
c-factor in the mBJ method. To this end, we reformulate 
the mBJ exchange potential [7] in an appropriate form 
as follows: 



/ I0k(r) 



, (1) 



where c, v^f 7 , t a , and p a are the c-factor, multiplicative 
Becke and Johnson potential [9], kinetic energy density, 
and electron charge density, respectively. The c-factor 
linearly depends on the square root of the average of 



Pa 



as follows [7]: 



1.023 




dr' -0.012, 



(2) 



where the two 1.023 \Jbohr and -0.012 (dimensionless) 
values were fitted by Tran and Blaha [7] in order to re- 
produce the experimental bandgaps of several com- 
pounds. From Equation 1, it is clear that if c = 1, then we 



have, v™^(r) = v® G (r) which is consistent with that of 
[7]. The c-factor in many cases can be well reproduced 
within a self-consistent calculation by the self- 
consistently converged charge density, p a . However, for 
the y-Al 2 0 3 compound, the self-consistently reproduced 
c-factor is not calculated to be so large to be used for 
predicting the bandgap close to the experimental value. 
The c-factor may not be large enough for our case be- 
cause the electron charge density for the defective spinel 
□ 2 |A/ 21 i0 32 around the Al vacancy, as indicated by the 
empty box (□), is very small. Therefore, due to the va- 
cancy, the c-factor which linearly depends on the square 

root of the average of ^-y^ is self-consistently calculated 
by Equation 2 to be smaller than an actual and necessary 
c value for reproducing the bandgap of the y-Al 2 0 3 . 
Hence, we have increased the c-factor to overcome the 
low electron charge density that originated from the va- 
cancy for our case. It is worth to note that, in practice, 
the c-factor cannot be unlimitedly increased. Thus, the 
experimental bandgap cannot be always exactly repro- 
duced by increasing the c-factor in the current version 
of the mBJ method for every case without considering 
the electron charge density variation and bonding nature 
of the case. Indeed, there can be a critical value for the 
c-factor, as reported in [10]. For larger c values than the 
critical c value, the bandgap may be drastically decreased 
to a meaningless value far from the experiment. This 
shows that the mBJ potential may result in an incorrect 
electron charge density, if one uses a very large c-factor. 
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The latter point was perfectly demonstrated [10] by cal- 
culating the electric field gradient, as an extremely sensi- 
tive physical quantity to the shape of the valence 
electron charge density, versus the c-factor for the Cu 2 0 
compound. Consequently, the c-factor, as a measure for 
the electron charge density, should be carefully exam- 
ined and optimized at least for those cases having low 
electron charge density. The low charge density may ori- 
ginate from some hollow spaces due to the vacancies, as 
what existed in the y-Al 2 0 3 case, or from weak van der 
Waals bonds, as what existed in the fcc-C 60 fullerite. 
Therefore, from the above discussion, we would, as a 
corollary result, anticipate that the c-factor may not be 
also well produced by the regular mBJ method for the 
case of fcc-C 60 . Thereby, one needs to first optimize the 
c-factor for this case as well. The latter anticipation 
needs more elaboration which is out of the scope of this 
work. 

Charge density 

A graphical description of the spatial features of the 
electronic structure has been presented by means of the 
electronic density [14]. Figure 4 presents the electronic 
charge density graph for the defective spinel structure of 
y-Al 2 0 3 which is calculated using the mBJ exchange po- 
tential. Clearly, the electronic cloud is mostly distributed 
around the oxygen nuclei. Moreover, in the Al-O inter- 
nuclear line, the electronic cloud is polarized toward the 
aluminum atoms. Neither the Al-Al nor O-O inter-nu- 
clear lines show a polarized electronic cloud. This indi- 
cates that a highly polar bond is formed in the Al-O 
line, in y-Al 2 0 3 , as expected from a large electronegativ- 
ity difference between Al and O atoms. These results are 
consistent with previous reports (e.g., [6]). It should be 
pointed out that electron depletion is observed at the 
aluminum vacancy sites in Figure 4. Observation of such 
an electron depleted region is expected because a stoi- 
chiometric y-Al 2 0 3 structure itself is an electrically neu- 
tral solid material, and inclusion of an ionized aluminum 
vacancy (e.g., Vai" point defects) makes the electric 



charges unbalanced. Moreover, the stoichiometric 
y-Al 2 0 3 structure is a perfect insulator, and formation 
of charge-free carriers is not anticipated from our calcu- 
lation. In addition, it is observed in Figure 4 that distri- 
bution of the valence charge density is not uniform 
around different Al lattice sites. These results do not 
completely confirm the previous results that distribution 
of the valence charge density is uniform around the Al 
lattice sites, regardless if they are occupied or vacant 
[14]. We believe that our findings may be a little bit in 
better agreement with the y-Al 2 0 3 material because 
introduction of a vacancy affects the aluminum- oxygen 
interactions and therefore slightly disturbs the charge 
density distribution. 

Figure 4 qualitatively shows that the Al atoms located 
at the octahedral positions, in contrast with the oxygen 
atoms, are covalently bonded to the O atoms. Al atoms 
are far from each other and thereby cannot be strongly 
bonded to each other. These observations are in agree- 
ment with the quantitative ionic charge analysis [14]. 
Therefore, the Al-O constitutes the most important 
bonding in the y-Al 2 0 3 crystal. This implies that if the 
electron charge density of the y-Al 2 0 3 crystal is sub- 
tracted from the superposition of its atomic densities, 
the Al-Al and O-O charge densities can almost cancel 
each other. This is consistent with the difference be- 
tween the electron charge density of the y-Al 2 0 3 crystal 
and the superposition of its atomic densities [14]. There 
are two Al vacancies in the supercell, as stated in the 
'Crystal structure' subsection and as shown in Figure 1, 
and two positions with two different point groups for 
placing them, octahedral (O) and tetrahedral (T) posi- 
tions. Either both of the Al vacancies can be positioned 
in the octahedral (tetrahedral) and OO (TT) sites or one 
in the O site and the other in the T site, OT, as consid- 
ered in [5]. Depending on the vacancies positions, 14 
nonequivalent, 4 OO, 6 TO, and 4 TT configurations 
can be constructed for the system [5]. Total energy 
DFT calculations were performed using the ultrasoft 
Vanderbilt pseudopotential as embodied in the VASP 
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code for these 14 configurations [5]. The minimum en- 
ergy was found for an OO configuration, followed by an 
OT, and then a TT configuration with energy differences 
of 0.160 eV/(Al 2 0 3 ) and 0.137 eV/(Al 2 0 3 ) [5]. These 
small energy differences show that the total energy is 
not completely indifferent to the Al positions. This veri- 
fies that the most crystal bonding energies originate 
from Al-Al bonds. However, by neglecting these small 
energy differences, we would state that the total energy 
is almost insensitive to the Al positions, which is what 
would explain the disorder in the Al sublattice in agree- 
ment with [14]. Equation 1 can be equivalently rewritten 
in terms of Becke-Roussel (BR) exchange potential [15] 
as follows [7,10]: 

The BR exchange potential tries to model the Cou- 
lomb potential created by the exchange hole and 
approaches asymptotically to the exact Kohn-Sham ex- 
change potential [10]. The mBJ potential as given in 
Equation (3) consists of two terms - the BR term and 
the \ftfp term [10]. In addition to the total electron 



charge density, we have calculated the electron charge 
density generated only by the BR term, as shown in 
Figure 5a. The BR term is calculated by setting c = 2/3 in 
Equation 3. Furthermore, we have calculated the elec- 
tron charge density generated only by the second term, 
yjtjp term, as shown in Figure 5b. The y/t/p term is 
calculated by setting c = 0 in Equation 3. The BR charge 
density, as shown in Figure 5a, is almost spherical, while 
the charge density is aspherical for the y/t/p term, as 
can be seen in Figure 5b. This result is in agreement 
with [10]. Exchange potentials that corresponded to the 
BR and y/ t/p terms are similarly calculated. The vj^(r), 

y^, and (r) exchange potentials are given in two- 
dimensional plots, as shown in Figure 6a,b,c, and in 3D 
plots, as shown in Figure 7a,b,c, respectively. The results 

show that the v^(r) is attractive (negative), while y^^y 
is repulsive (positive) which is in complete accord with 
[10]. The total exchange potential, v^(r), (Figure 6c) 
can be obtained by applying Equation 3 on the vJJ^(r) 

(Figure 6a), and ^Jy^ (Figure 6b) contributions for our 

optimized o factor, i.e., c = 1.8. For example, for the red 
regions of Figure 6a,b,c, we can verify the latter point, 
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+0.0699 
+0.3029 
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+5.6900 
•■24.6604 



(a) 



| Scale: A n(r) | 




(b) 

Figure 5 BR and \j*~[ P terms corresponded to the electron charge density of y-AI 2 0 3 . (a) BR and (b) 

electron charge density, as given in Figure 4, indicating the position of the aluminum vacancy for y-AI 2 0 3 . 
by setting c = §, and \ r / P is obtained from Equation 3 by setting c = 0. 



/ p terms corresponded to the total 
term is obtained from Equation 3 
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Figure 6 Contour plots for v B x R a {r) and y^P) contributions of the total v™ B J(r) exchange potential of y-AI 2 0 3 . Two-dimensional plots for 
(a) v^ 0 (7) and (b) contributions of (c) the total v™$ J (7) exchange potential for y-AI 2 0 3 . v b x r 0 (j) is obtained from Equation 3 by setting c = \ , 
and J r / P is obtained from Equation 3 by setting c = 0. v^ J (7) is obtained from Equation 3 using the optimized c-factor, i.e., c= 1.8. 



viz., 213.5565=1.8 x (-260.4229) + (3 x 1.8 - 2) x \ x 

y/§ x 256.8886 = -214.9687. Therefore, it is now clear 
that the bandgap should be increased by increasing the c- 
factor: the signs of the BR and \ftfp terms are opposite, 
and BR is weighted by c, whereas yjtjp is weighted by 3c 
[10]. The repulsive potential, yjtjp is proportional to 
\fljp. Therefore, the \ftjp term in the mBJ potential will 
be of significant importance if the electron charge density 
is low. For our case, charge density is low in the regions 
near the vacancies. Thus, contribution of the repulsive 
term is important for the A1 2 0 3 case, as can be seen in 
Figure 6b, and should not be neglected. 

Conclusions 

We have successfully simulated the y-Al 2 0 3 compound 
based on the DFT method using the full potential aug- 
mented plan waves plus local orbital method, as em- 
bodied in the WIEN2k code, and applied the mBJ 
exchange potential on this system to predict its bandgap 



more precisely. We found that the Al-O has a highly 
polar bond in this compound which is consistent with 
previous reports. We showed that distribution of the 
valence charge density is not uniform around different 
Al lattice sites. Besides, the calculated partial DOSs indi- 
cated that majority of the valence electronic charges cor- 
respond to the p orbitals of oxygen atoms which are 
consistent with those of previous reports. However, we 
showed that different oxygen atoms have different con- 
tributions in the valence electronic charges. Contribu- 
tions of oxygen atoms which are closer to the vacancies 
in y-Al 2 0 3 dominate. A charge-free region is observed 
at the aluminum vacancy site which confirms the charge 
neutrality and insulating behavior of the stoichiometric 
y-Al 2 0 3 structure. Our result shows that mBJ can sig- 
nificantly improve the electronic structure of the system 
if a suitable c-factor is used. A direct bandgap of 8.02 eV, 
which is very close to the experimentally measured value 
of 8.7 eV, was obtained at the T point by adjusting the c- 
factor internally to a value of 1.8. This c-factor value can 
be utilized for the correct estimation of the electronic 
and optical properties of y-Al 2 0 3 compound based on 
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(a) 




(b) 




Figure 7 Three-dimensional plots for vj* (r) and ^J^j^ 
contributions of the total v^ J (7) exchange potential of 
Y-Al 2 0 3 . Three-dimensional plots for (a) v B x R 0 (j) and (b) 

contributions of (c) the total v^ J (T) exchange potential for y-AI 2 0 3 . 
v B x R 0 {f) is obtained from Equation 3 by setting c = |, and yj*~/p is 
obtained from Equation 3 by setting c = 0. v^ J (7) is obtained from 
Equation 3 using the optimized c-factor, i.e., c= 1.8. 



the full potential augmented plan waves plus local or- 
bital DFT method using the mBJ exchange potential. 
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